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Rate equations are used to model spectroscopic observation of incoherent energy transfer in light- 
harvesting antenna systems based upon known structures. A two-parameter two-dimensional model 
is proposed. The transfer rates obtained, by matching the fluorescent decay, are self-consistent 
within our model. 
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I. INTRODUCTION 



We have had a reasonably complete picture of the bac- 
terial light-harvesting (LH) system recently P10]. Both 
the inner antenna, LH1, and the outer antenna, LH2, 
are assembled from the same modules to form rings. 
Each module consists of two short a-helical polypep- 
tides coordinate one carotenoid and three bacteriochloro- 
phylls (BChls). The LH2 is composed of 9 units, for 
Rhodopseudomonas acidophila (^], resemble a cylinder, 
with an inner diameter 36A and an outer diameter 68A, 
while the LH1 is composed of 16 units, for Rhodospir- 
illum rubrum 0, in order to accommodate the reaction 
center (RC). The later has an outer diameter 116A and 
a central diameter 68A However, the exact numbers of 
both complexes are variable [QjHJ^]. 

Furthermore, the LH2 B850 BChl a form a complete 
over-lapping ring in a hydrophobic environment, which 
reduces the dielectric constant, while the B800 BChl a 
are well separated and are in a polar environment. When 
a BChl molecule is excited by light, the energy can reach 
equilibrium within about lOps Q. A LH2 can function 
as a storage ring to store the excited singlet state energy 
for about llOOps. However, the energy will transfer to 
other rings before decaying. The hopping of energy con- 
tinues from one ring to another one until a LH1, which 
contains the RC, is finally reached. The total trip lasts 
for about 5 to 50ps Apparently, there is a com- 

petition between energy relaxation and energy transfer. 



Historically, relatively few physicists have tackled 
problems of photosynthesis. Notably, Montroll used ran- 
dom walk concepts to model energy transfer amongst an- 
tenna rings on a lattice by considering its first passage 
time [jio). Later, Hemenger et al. proposed a more real- 
istic model by taking inhomogeneous transfer rates and 
trapping of RCs into account jnj. Interestingly, it is 
Pearlstein's work which is most often cited in the liter- 
ature [ p^[ . In the mean time, almost all experimental- 
ists try to find some explanations for their spectral data. 
However, due to lack of precise geometrical information 
most efforts are in vain. 

Progresses in physics are often made along the line 
structures - energy - dynamics. A goal of researches 
nowadays is to find the relation between structural and 
spectral information obtained, expecting that the func- 
tion of photosynthesis will be explained in terms of its 
structure, and further drawing inferences from the model 
by applying methods of mathematical or numerical anal- 
ysis. Recently Timpmann et al. used a rate equation 
model to describe energy trapping and detrapping by the 
RC p~3| ] . However, their antenna has no structure. Skala 
et al. also carried out a series of investigation by analyz- 
ing the spectrum of a more r ealistic LH1 model 
However, their model is incompatible with the recent 
structural finding. In this paper we established a two- 
parameter model based on recent structural data. 
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II. MODEL 



With the known periodical structure, shown in Fig.[l], we can built, from chemical rate equation, the following 
phcnomenological model of energy transfer, 



^ = k Ai - (k" + k E )E 



dAi 
~dT 



kA 16 - 2kA x + kA 2 - k A x + k E , 



(1) 
(2) 
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dt 
dA 16 

dt 



feA„_i - 2k A n + kA n+1 , 



= kA 15 - 2kA w + kAx 



2. .Ah , 



(3) 
(4) 



in which A n s denote the excited BChl dimer, E = P*BH is the excited state, with B representing the chlorophyll 
monomer within the RC, and P* is the excited special pair of BChl molecules. It is a set of 17 coupled linear 
differential equations. The symmetry of this system is broken due to k ^ k . A similar model has been proposed 
by Skala et al. [ fl5| |. However, the RC and the antenna ring are connected only at one site in the present model, 
corresponding to the recent experimental observation. 
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FIG. 1. Schematic plot of LH 1 and definition of symbols used. 



In the homogeneous case with the same transition rate amongst the units, the characteristic polynomial of the 
above rate-constant-matrix can always be expressed as 



p pi p2 p3 p4 

^16 — - r 16- r 16- r 16- r 16 



(5) 
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with 



f* = s + 2k , (6) 

pfg = s 2 + Aks + 2k 2 , (7) 

pfg = s 4 + 8fcs 3 + 20fc 2 s 2 + 16fc 3 s + 2fc 4 , (8) 

P 4 g = s 10 + fa + k " +k ' + 18fc)s 9 + 



(k'k E + 18kk E + 18kk" + Wkk + 134fc 2 )s 8 + 
2(8k'k E + Q7kk E + 67kk" + h2kk + 266fc 2 )fcs 7 + 
2(52k'k E + 266fcfc B + 266fc/c" + 176fcfc' + 605fc 2 )fc 2 s 6 + 
2(176k'k E + 605fcfc B + 605fc/c" + 330fc/c' + 786fc 2 )fc 3 s 5 + 
12(55k'k E + 131kk E + 131fcfc" + 56kk + 91fc 2 )fc 4 s 4 + 
4(168k'k E + 273kk E + 273kk" + 84kk + 86fc 2 )fc 5 s 3 + 
8(42k'k E + AZkk E + 43fcfc" + 8kk + 4fc 2 )fc 6 s 2 + 

2(32k'k E + \Qkk E + Wkk" + kk')k 7 s + 2k 8 k'k E , (9) 

which is a consequence of the master equation used, and is independent of the detail geometrical symmetry. The 
mode controlling the decay to the RC is within P 4 6 , since P^ 6 , P 2 6 , P 3 6 do not contain k , k and k E . However, 
all four parts will be influenced by the change of k. If one solves this set of differential equations by applying the 
Laplace transformation method, one finds the solution divides into four distinct groups of decay channels, namely, 
A5-A13; E-Ai-Ag; A3-A7-A11-A15; A2-A4-A6-A8-A10-A12-A14-A16. Because the matrix of rate constants is hermitian, 
all eigenvalues are negative. Furthermore, no eigenvalues are degenerated, in contrast to Skala's model which posses 
too high degree of symmetry pi . Letting k = k does not results in additional factorizability although the symmetry 
of our model is restored. At k = k =0, Pie becomes 

s(s + 2k) 2 {s + 4fc)(s + k E ){s 2 + Aks + 2k 2 ) 2 (s i + 8ks 3 + 20fcV + 16fc 3 s + 2fc 4 ) 2 . (10) 

It contains a zero eigenvalue, which signals the existence of a steady-state solution, as should be happened without 
the decay to the RC. Degeneracy of eigenvalues is introduced as the transition to the RC is decreased. 



III. SPECTROMETRY COMPARISON 



We can verify our model against experiments: The 
pump-probe spectroscopy measures the difference be- 
tween two beams, with 

AD = Ae A J2A n + Ae E E, (11) 

n 

being the signal measured. The symbol Aes are the 
differences in dielectrical constants between pump and 
probe beams of the corresponding pigments. By choosing 
the pump and probe laser frequencies, we can selectively 
detect the population changes of ^2 A n or E. Summing 
over Eq. (|l|)-(Q) we know that the decay of the total 
population should be d(^2 A n )/dt = —k A% + k E. The 
measured charge separation rate is k E « 3.57 x 10 11 s~ 1 
at room temperature, and increases by 2 to 4 times from 
300!^ to 10K depending on the species chosen ]l7| , |l8| |. 
The ratio of the forward and backward transition to the 
RC is know to be about 25% for an open RC, i.e., 
the RC BChl dimer (P) is reduced and the iron quinone 
electron acceptor is oxidized; 40% for pre-reduced RC. 
The back-trapping rate can, in principle, be estimated 



from k jk — exp(— AG/k E T), with AG the free-energy 
gap between A\ and E is estimated from their absorp- 
tion peaks, k E is the Boltzmann constant, and T is the 
absolute temperature. However, the measured absorp- 
tion peaks of the excited RC are broad and imprecise [|| . 
We do not know the absolute values of k or k experi- 
mentally since it is difficult to tune the laser frequency 
to distinguish A n from E. Nor do we know the tran- 
sition rate between A n s because transition between the 
same species cannot be measured directly. Furthermore, 
at room temperature, energy equilibrium within the an- 
tenna interferes with the trapping process. Therefore we 
have taken k and k as parameters and fit the slow mode 
of fluorescence decay of excited population observed, i.e. 
200ps H^Q- Thus, the absolute value of the largest 
eigenvalue should be about 3/200ps = 1.5 x 10 10 s _1 . 
A computer code is written to scan all combinations of 
k and k , with k — k /5, for the largest eigenvalue 
to be smaller than —1.5 x 10 10 s _1 between — 10 8 s~ 1 to 
— 10 15 s _1 . Interestingly, we find all possibility happened 
at k — k and k > 6.97 x 10 11 s _1 for k = k /5. Presum- 
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ably, it is an extremum of Pf 6 . At the lowest k, we can 
match the required 200ps decay whose curve is plotted 
at Fig.§ If k" = fc'/4, we obtained k = 7.25 x lO 1 ^" 1 . 
That k has to be equal to k might sound peculiar in 
view of the geometrical distance between Ai and RC is 
less than the distance between RC and other A n s [Q. 
However, the species for donor and acceptors are differ- 
ent at these two cases. There are possibilities that the 
final hopping rate are still the same. 

signal Fig. 2 Ting 
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2. Numerically calculated pump-probe signal from 



at k = 6.97 x 10 11 , k E = 3.57 x 10 11 , k = k. 



FIG 
Eq. (0) 

k — k /5, Aea = 1, Ae_B = 0. The initial condition is 
A 5 = 0.2, A 7 = 0.4, Aio = 0.3, A 12 = 0.1, while other sites 
are not excited at t = 0. 



The transfer of excitation energy requires coupling 
between the emitting molecule and the ground state 
molecule. At an intermolecular separation involved be- 
tween 10 A to 10CL4, long-range resonance transfer of elec- 
tronic excitation arises from coupling between the tran- 
sition dipoles of the donor and the acceptor, which is the 
Foster theory j^j2(J. Since the BChl Q y dipoles lie in 
the same plane, we have 



k(R) cc — ( 

Tp R 



(12) 



in which i?n, measures transfer efficiency, is the Foster 
radius, van Grondelle gave Rq — 9QA for the BChl 
875 to BChl 875 energy transfer and a fluorescence life 
time, Tp, about 3000ps or slightly higher pl| , p2| . If a 
putative separation distance between interacting BChl 
a dimers ~ 17. 5A is used Q we obtain an estimation 
of k w 6.17 x 10 s . This number is about an order 
of magnitude higher than the value obtained from our 
model. However, the pairwise energy transfer is about 
lps according to our calculation ||. On the other hand, 
from the value of k obtained here, by fitting the 200ps de- 
cay as well as the tf, we estimated the Foster radius to be 
26. 8A. This result is consistent within our model since 
we assume only nearest neighbour transition. Further, 
since we put the population at the antenna at t = for 
our calculation, the rising time is infinitely short, instead 
of having some instrumental limits as observed experi- 
mentally. Although the light wave length is much larger 



than the ring size, the ring still might receive energy in 
localized form by energy transfer from other rings as the 
initial condition we used in Fig.||. Table I provides a 
list of all eigenvalues and corresponding amplitudes ob- 
tained from our model. From the table, we found that 
the largest eigenvalue mode is important, not only for its 
large separation from the other eigenvalues but also for 
its corresponding large amplitude. 



amplitude 


eigenvalue 


-0.0070 


-2.9707 


x 10 12 


-0.0006 


-2.7443 


x 10 12 


0.0000 


-2.6819 


x 10 12 


0.0221 


-2.4614 


x 10 12 


0.0000 


-2.3797 


x 10 12 


-0.0200 


-2.0143 


x 10 12 


0.0000 


-1.9275 


x 10 12 


0.0124 


-1.4844 


x 10 12 


0.0000 


-1.3940 


x 10 12 


-0.0086 


-9.5864 


x 10 11 


0.0000 


-8.6054 


x 10 11 


-0.0159 


-5.5606 


x 10 11 


0.0000 


-4.0829 


x 10 11 


0.0034 


-3.8570 


x 10 11 


-0.0716 


-1.4875 


x 10 11 


0.0000 


-1.0611 


x 10 11 


1.0858 


-1.5107 


x 10 10 



TABLE I. Eigenvalues 
time-resolved experiment 
respond to the lifetimes 
fluorescence or bleaching 



for LH 1 for Fig.|2[ In a 

the relaxation rates cor- 
observed from antenna 
kinetics. 
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We have also introduced inhomogeneity, due to geo- 
metrical distortion, into the rate constant. However, even 
at large distortion, the basic character of the spectrum is 
not altered considerably. If the criteria for k = k can be 
established, we can further reduce the free parameters in 
our model. 



IV. CONCLUSION 

In summary, a physicist's approach p3| of incoherent 
energy transfer within an antenna ring is taken by con- 
sidering a two-parameter two-dimensional model. This 
model differs from the one presented by Skala et at. The 
reality might be somewhere between these two models. In 
our model, we numerically found k has to be equal to k . 
Furthermore, we are able to calculate analytically some 
of the eigenvalues and demonstrate explicitly that there 
is a mode for decaying to the RC and other three modes. 
However, this result of mode separation depends upon 
the exact number of unit involved in the ring. Therefore 
should not be important. Perhaps we should interpret 
the finding as: P-j^, P^, Pf 6 are redundant, since Pf 6 
contains ks which should be important. A ring of 16 
units is huge. The only purpose for such a large antenna 
is to accommodate the RC. 

Finally we remark that it is possible to extend a two- 
dimensional random walk model of energy transfer into a 
quasi-three-dimensional one, in view of the recent struc- 
tural finding, with a recent result of random walk on bun- 
dled structures by Cassi and Regina pij . Furthermore, 
this theoretical result should be able to be verified exper- 
imentally using its spectral dimension by measurements 
involving diffusion processes such as time-resolved spec- 
troscopy of nearest-neighbours energy transfer. Other 
light-harvesting models and mechanisms are under fur- 
ther investigation. 
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